Gravitino Dark Matter with Stop as the NLSP 
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Abstract. Gravitino is one feasible candidate for the dark matter in supergravity models. With its 
couplings being suppressed by the Plank mass, gravitino interacts very weakly with other particles, 
making its direct detection, or production at colliders practically impossible. The signatures of this 
scenario mainly arise from the next lightest supersymmetric particle (NLSP) which is metastable 
due to the very weak coupling. There are many possibilities for the NLSP, and here I will review 
the feasibility of stop particle as the NLSP and discuss its phenomenology. 



PACS. 12.60.Jv Supersymmetric models - 95.35.+d Dark matter 



1 Introduction 

Even with many astrophysical evidence pointing to the 
existence of dark matter, we have not yet known what 
the dark matter is. On the other hand, studies in par- 
ticle physics suggest that it could be a new undiscov- 
ered particle from beyond the standard model that is 
massive, neutral and stable (or very stable with life- 
time longer than the age of the Universe) . In addition, 
because of the so far nil results of direct detection 
searches, the dark matter particle can only interact 
weakly (or very weakly) with ordinary particles. 

In supersymmetric models with i?-parity conserva- 
tion the lightest supersymmetric particle (LSP) is sta- 
ble. If the LSP is neutral it could be a candidate for 
dark matter. Indeed in the minimal supersymmetric 
standard model (MSSM), for example, we have neu- 
tralino and sneutrino as candidates for dark matter. 
There have been lots of studies on these possibilities, 
and it would be impossible to list all of the papers on 
these subjects here. While 'left handed' sneutrino dark 
matter is already excluded, neutralino is still a feasi- 
ble and appealing candidate for dark matter. Includ- 
ing gravity in supergravity models, gravitino is another 
neutral particle that has been shown as a suitable can- 
didate for the dark matter, see e.g. [2j[3]. The strength 
of gravitino interaction with other particles is sup- 
pressed by the Planck scale and as a consequence the 
next lightest suppersymmctric particle (NLSP) would 
be metastable, assuming that i?-parity is conserved. 
The phenomenology of this scenario depends much on 
what the NLSP is. Depending on the supersymmetric 
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model, it could be any of the supersymmetric particles 
beside the gravitino, each with its own phenomenol- 
ogy. Here wc would like to discuss the scenario with 
the lighter stop particle as the NLSP p] . We will show 
that it has distinct phenomenology, and discuss its im- 
pact on cosmology and supersymmetry searches at col- 
liders. 



2 Properties of Stop NLSP 

Let us start by looking at the stop mass matrix. In the 
left-right basis the mass matrix has the following form 
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Diagonalising the matrix, we get the mass eigenvalues 
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Wc see a seesaw mechanism here, i.e. one mass is sup- 
pressed while the other is enhanced by the off-diagonal 
terms in the mass matrix. Thus, if the off diagonal 
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terms are large, which is possible because mt is large, 
we could get a light stop. However, to get a stop which 
is lighter than other supersymmetric particles, in par- 
ticular in models with universal soft masses, it is not 
enough to have large mt but we need Aq to be also 
large. 

The light stop t\ would eventually decay to grav- 
itino G, which is assumed to be the LSP. Provided that 
it is kincmatically allowed, the dominant decay mode 
is the 2-body decay t\ — > G + t 
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If the mass gap A r , 



— m~ is less than m t but 
still greater than mw + ^6, the 3-body decay mode 
t\ — > G + W + b becomes the dominant one. It consists 
of three channels: 




The 3-body decay rate can be approximated by the 
following phase space expression 
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Note however that for the numerical calculation, we 
use a more complete analytic formula. The numerical 
results for the stop lifetime are shown in Fig. [TJ The 
first panel shows the lifetime r when 2-body decay is 
dominant. We see that the lifetime can be up to ~ 
10 9 s. Note that r decreases as A m increases. The fig- 
ure only shows results with m~ > 10 GeV. For smaller 
gravitino mass we could get lifetime shorter than 10 3 s. 
The second panel shows the lifetime through 3-body 
decay, and we see that the stop can live up to ~ 10 15 s. 

Of course it is possible that the mass gap between 
stop and gravitino is even smaller than mw, in which 
case the 3-body decay channels above are not available. 
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Fig. 1. Stop lifetime as a function of gravitino mass m~ 
with fixed values of m~ tl = 200, 300, 400, 500 GeV respec- 
tively from left to right. 



In this case we have to consider the 4-body decay tj — > 
G + b + (qq,lv). We can estimate the lifetime using 
phase space 



A-body « 10~ 30 GeV-V 1 ' 



\Amf ({Amf - mlf' 2 
(7) 

We shall not forget though, that stop can decay through 
loop processes as well, such as the following VF-loop 
diagram 




However, this process is suppressed by large , es- 
pecially in the models that we consider, and the small 
CKM mixing Vckm- So the bottom line is, depending 
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on A m and the model we use, the stop lifetime could 
be very long, even longer than the age of the Uni- 
verses (O(10 17 s)). However, the longer the lifetime the 
stronger the constraints, such as from the black body 
spectrum of CMB. For simplicity we consider lifetime 
in the range of 1 — 10 8 s, and assume that the cosmo- 
logical constraint is coming only from BBN. 

A long lived stop would hadronize. By taking ana- 
logue to heavy quark hadrons we can deduce the light- 
est hadronic states and approximate their lifetime. The 
lightest sbaryon is A~ = t%ud. There are also S~ + ' + ' = 

t\ (uu, ud, dd) sbaryons which would decay strongly into 
the lightest sbaryons, and the next ones are S~' = 
t\s{u, d) which decay semileptonically with lifetime r Ss 
10 -2 s. There are mesinos, with the lightest one is T° = 
t\u. We also have T + = t\d with lifetime r ~ 1.2 s, 
and T s = i\s with lifetime r ~ 2 x 10~ 6 s. Similarly 
the antistop would hadronize into the corresponding 
antisbaryons and antimesinos. So, in about one sec- 
ond after the big bang, only the lightest mesinos T°(*) 
and lightest sbaryons A— remain. 

3 Cosmological Implication 

It is well known that a metastable particle can alter 
the light elements abundances from the standard BBN 
prediction if it decays during or after the BBN time. 
The effects can be classified as follows: 

— Photodissociation: Electromagnetic showers from 
the decay can destroy light elements already formed 
by BBN. The products of this dissociation would 
then involve in the continuing nucleosynthesis pro- 
cesses. 

— Hadronic showers: There could also be effect from 
quarks and gluons produced by the decay. This ef- 
fect is of two types: 

— If the decay happened early, the hadronic show- 
ers could change n/p ratio, and we call it hadron 
injection effect. 

— If the decay happend around the beginning of 
the BBN, the hadronic showers could also dis- 
sociate light elements, especially obgj an d we 
call this hadrodissociation. 

The hadronic showers become relatively unimpor- 
tant if the decay happens after around 10 6 s. 

— Catalytic bound state effect: As pointed out quite 
recently by Pospelov [4] there could be another im- 
portant effect if the metastable particle is nega- 
tively charged. In this case it can form bound state 
with nuclei, hence lowering the Coulomb barrier for 
certain nucleosynthesis proccscss and introducing 
photonlcss final states for radiative capture reac- 
tions with much enhanced cross sections. 

Now, looking at the stop case, there are two signif- 
icant facts influencing its BBN effect. First of all, as 
a strong interacting particle, stop naturally has small 
relic density due to the large annihilation cross section. 
In turn, the photodissociation and hadronic showers 



effects would also be relatively small. Secondly, as al- 
ready discussed above, after hadronization only vi- 
and ?°W left. Because of the mass difference, T° and 
its antiparticle are more abundance than A~ by ~ 

O(10) during the early stage. The negatively charged 
A~ is further suppressed by: (1) pairing and subse- 
quent annihilation with y1~; and (2) quark exchange 

with ordinary hadrons (proton and neutron) into T . 
Note that A— would also annihilate with antibaryons 

although not as much as the AZ annihilation with 

T 

baryons. Thus at the time of BBN practically only the 
neutral mesinos around, and the bound state catalytic 
effect would become irrelevant. The BBN constraint 
has two parameters: the density of the metastable par- 
ticle before decay and the lifetime (with fixed baryon 
density). The stop lifetime depends also on the grav- 
itino mass, which we assume to be a free parameter. 
As long as the lifetime is not too long (i.e. ^ 10 s s) we 
should be able to satisfy the BBN constraint. 



4 Collider Signatures 

We can also look at the production and detection of 
metastable stop at colliders. Provided there is enough 
energy, stop and antistop could be pair produced, and 
(assuming metastable) they would hadronize before 
passing the detector. 

There would be neutral as well as charged shadrons 
(both sbaryons and mesinos). As the stop shadrons 
propagate they could undergo charge exchange as well 
as baryon exchange with the background nucleons, for 
example T + (p,n) — > (Af,Sf) + nir. We estimate 
that only about 1/16 of the produced stop-antistop 
pairs yield clear signal of opposite charged pair. 

The charged sbaryon would appear as a charged 
heavy particle in the muon detector. Looking for 'slow 
muon' signal and combine with stop production cross 
section, one can set the metastable stop mass lower 
limit. There has been analyses by CDF on the Teva- 
tron Run II data in order to set the mass limit. The 
most recent limit is 250 GeV [5], although at the time 
we did our analysis the limit was mj > 220 GeV [6]. 
Both numbers are not published yet, and the limit 
from the Tevatron Run II could still be changed. 



5 Stop NLSP in Specific Models 

The next obvious question is: can we realize this stop 
NLSP scenario in a specific supcrsymmetric model? In 
the CMSSM it turns out to be not possible. The reason 
is because three constraints: stop being the NLSP, the 
Higgs mass lower bound, and the stop mass lower limit; 
do not have an overlap allowed region. 

In the NUHM model [7j, which has two more free 
parameters compared to the CMSSM, we have a better 
luck and found a small region that is still allowed for 
mi > 220 GeV. It is shown in Fig. [2j Note however 
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m 1/2 = 600, m = 500, m A = 1400, A = 2100, tan (3 = 10 
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Fig. 2. NUHM. 

that this region barely exists if we take 250 GeV as 
the stop lower mass limit. 

This does not mean that this scenario is in general 
almost completely excluded. Let us take a look at com- 
pressed supersymmetry model [8], for example. This 
model is motivated by SUSY little hierarchy problem, 
i.e. fine tuning between |/i| 2 and m 2 H which is required 
by electroweak symmetry breaking condition. This fine 
tuning could be aleviated if n%H u is not large. As we 
know from the RGE, the biggest contribution to m 2 H 
comes from Ms. Thus, this suggests a model with small 
gluino mass. Now, light gluino also leads to light stop, 
as can be seen through the RGE. In [8] , the author still 
considered neutralino LSP scenario without gravitino. 
However, if we choose a bit smaller M3 or a larger but 
still modest we would be able to get a stop lighter 
than neutralino. Note however, that the Higgs mass 
is strongly correlated to the stop sector, so a detail 
study would be required to reach conclusive statement 
on this model. 

Perhaps an easier way to satisfy all the constraints 
is to lift up the neutralino mass instead, i.e. by choos- 
ing large bino mass M\ in a nonuniversal gaugino masses 
model. With large neutralino mass, it would be easy 
to find a set of parameters with stop lighter than neu- 
tralino but still relatively heavy. 



6 Metastable Neutralino 

Looking at Fig. [2J we notice that x? — > i\ + t is not 
kinematically allowed, hence we could have a long lived 
neutralino as well. This has implication on cosmology 
as well as in collider physics. 

In the early universe, neutralino and stop coexist 
for sometime, with relative density fi x ~ ^ J?^ . The 
neutralino could decay directly to gravitino, or to stop 
which then decay to gravitino - a cascade decay. All 
of these decays would affect the BBN, and therefore 
would complicate the analyses. 



Table 1. Supersymmetric mass spectrum for an NUHM 
model with allowed stop NLSP. 
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At colliders, there would be both missing energy 
signal of neutralino and slow muon signal of stop shadron. 
This suggests another benchmark point study for the 
LHC, with typical mass spectrum shown in Table [TJ 

7 Concluding Remarks 

In this review, we have shown that stop NLSP with 
gravitino dark matter scenario is very interesting from 
phenomenological point of view. Stop naturally has 
low relic density due to its strong interaction, making 
it possible to satisfy the BBN constraint. Metastable 
stop would hadronize, and the heavier shadron states 
would decay to the lightest ones. At the time of BBN 
practically only the lightest neutral mesino left, hence 
there is no EM bound state effect. 

This scenario is not feasible in the CMSSM because 
of the combined constraints from the stop NLSP, the 
stop mass and the Higgs mass lower bounds. In the 
NUHM this scenario is still barely possible. Nonuni- 
versal gaugino models with light gluino and/or heavy 
bino should revive the feasibility. 
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